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Abstract 

The difference between the structures of jets produced in heavy-ion and hadronic 
cohisions can best be exhibited in the correlations between particles within those jets. 
We study the dihadron correlations in jets in the framework of parton recombination. 
Two types of triggers, vr^ and proton, are considered. It is shown that the recombi- 
nation of thermal and shower partons makes the most important contribution to the 
spectra of the associated particles at intermediate pT- In pp collisions the only sig- 
nificant contribution arises from shower-shower recombination, which is negligible in 
heavy-ion collisions. Moments of the associated-particle distributions are calculated to 
provide simple summary of the jet structures for easy comparison with experiments. 
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1 Introduction 

Recent experiments at RHIC have revealed extensive information on the effects of the dense 
medium on hadron production at large transverse momentum (pt) [1, 2, 3]. The suppression 
of back-to-back correlation relative to the same-side correlation is a strong indication of 
substantial energy loss suffered by hard partons propagating through the medium [4]. It 
implies then that the jets detected in heavy-ion collisions are produced mainly near the 
surface of the medium so that the hard partons are less attenuated by the jet quenching 
effects [5]. If that is indeed the case, then a simple fragmentation model would predict that 
the structure of jets produced in heavy- ion collisions should be basically the same as that 
of jets in pp collisions. That similarity has been shown to be absent in the data of more 
recent experiments [6]. The aim of this paper is to investigate the jet structure by examining 
the two-particle correlation within a jet produced by nuclear collisions at high energy. The 
difference from pp jets is naturally caused by the presence of thermal partons in the jet 
environment, which we shall take into account. 

The framework in which we shall study hadron production at large pr is parton recom- 
bination, which has been shown to explain some features of the data where fragmentation 
fails [8]- [11]. The anomaly associated with species dependence of the Cronin effect in d-l-Au 
collisions has also been resolved in the recombination picture [12, 13]. The main compo- 
nent that is new in this series of work is the inclusion of shower partons generated by hard 
partons. Whereas the recombination of two shower partons forms a hadron that can be 
identified with the fragmentation product in the conventional approach, the recombination 
of shower partons with thermal partons yield hadrons in the intermediate pt region that are 
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totally new. Here we go a step further. To study dihadrons in a jet we must consider four 
partons that recombine to form two hadrons; some of those partons will be thermal in order 
to have enhanced yield. Dihadrons formed without thermal partons correspond to those 
found in jets in vacuum, and are suppressed compared to those hadron pairs that involve 
the participation of thermal partons. 

We shall consider two types of triggers, pion and proton, and calculate the distributions 
of the associated particles. Since we include the contributions from different species of hard 
partons, each having various flavors of shower partons, and since the four recombining partons 
have different momentum fractions that have to be permuted in the recombination formula, 
the combinatorial complication can result in a hundred terms or more. For that reason we 
limit our trigger to only tt"*" and p, which are sufficient to reveal the properties of the jet 
structures. Detailed comparison of our results to current data is, however, difficult, since 
the experimental trigger at this stage consists of all charged hadrons, as are the associated 
particles [4, 6, 7]. Separating triggers to mesons and baryons has not resulted in any pt 
distributions for the associated particles [14]. Nevertheless, some coarse comparisons with 
the data can be made, and our results will be shown to be reasonable. 

2 Single- and Two-particle Distributions 

The hadronization process that we consider in order to calculate the dihadron correlation 
in a jet is parton recombination. The formalism for single-particle inclusive distribution at 
high Pt is given in [11]. An essential ingredient in the hadronization process is the shower 
parton distributions (SPD), which give the probabilities of finding shower partons of various 
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flavors and momentum fractions in jets initiated by different hard partons [15]. Convoluting 
the SPD's, denoted by S^, with the hard-scattered parton distributions fi{k) in heavy-ion 
coUisions gives the shower component 

S{q,)=^J2 [dkkMk)Si{q,/k) , (1) 

i ■' 

where the sum is over aU hard partons, and ^ is the average suppression factor due to energy 
loss in a dense medium, found to be 0.07 for central Au-|-Au collisions at ^Js — 200 GeV 
[11]. For two shower partons in the same jet we have 

SS{q,, ^2) = e E / dkkm [Sl (I) , si (I) I , (2) 
where the curly brackets denote the symmetrization of the leading parton momentum fraction 



(3) 



.1-zJ '\1-Z2. 

This symmetrization is necessary, since either j or / shower parton may be the leading 
parton in the jet, and Eq. (3) is a way to ensure that momentum conservation Zi + Z2 < I 
is not violated. The SPD's are determined by solving the recombination formula for the 
fragmentation functions (FF) 



xD^ix) = / ^^{Siix,), S{{x2)}R^{x^,X2,x) . (4) 

J Xi X2 ^ ' 



where is the recombination function (RF) for the hadronization process j -\- j — > M. 

For the hadronization of a hard parton in vacuum it is unnecessary to consider the shower 
partons, since they recombine to form M, as in Eq. (4), to recover the FF, from which 
the SPD's are obtained. However, in the environment of thermal partons as in heavy-ion 
collisions the shower partons can recombine with the thermal partons, resulting in hadrons 



that are dominant in the intermediate px region because they benefit from the high density 
of the thermal partons as well as the higher momenta of the semi-hard shower partons. The 
invariant distribution of the thermal partons is parameterized by 

T{qi) =qi-T^ = Cq, exp(-gi/r), (5) 

where T is the inverse slope enhanced by flow. Since we have no model to describe the 
hydrodynamical evolution of the bulk medium at low transverse momentum, the parameters 
C and T are determined by fitting the low-pr data. They are found to be [11] 

C = 23.2 GeV-\ T = 0.317 GeV, (6) 

for central Au+Au collisions at midrapidity and ^/s = 200 GeV. 

With both the thermal T{qi) and shower distributions S{q2) known, it is possible to 
calculate the hadron distribution in the ID recombination formalism for the formation of a 

meson [16, 17] 

p dNM fdqidq2 , . , . , . 

p —r- = / Fqq'yqi,q2)RM{qi,q2,p) , (7) 

dp J qi q2 

where Fqq/{qi, is the joint distribution of a quark q at qi and an antiquark q' at q2- For 
pion production at high px, we put p in the transverse plane, abbreviate pt as p, ignore pion 
mass and write [11] 

where the RF [17, 18] 

i?.(?i,?2,p) = ^5f^ + ^-l), (9) 
p'^ \p p J 
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has been used. For heavy-ion coUisions Fqqi can be written in the form 

Fqq, ^TT + TS + SS , (10) 

where the possibihty of two shower partons from two different jets is ignored. It is the 
TS term in Eq. (10) that dominates at intermediate pr, while the TT and SS terms are 
dominant in the low and very high pt regions, respectively. 

Two-pion distribution can be obtained by a straightforward extension of the single- 
particle distribution given in Eq. (7), and one gets 



dN 



/ n — -^4(?l,?2,g3,?4)-R7ri(gi,?2,Pl)-R^2(g3,?4,P2), (H) 



PiP2dpidp2 pIpI J Vi=i * 
where sums over different combinations of partons contributing to the two RF's are not ex- 
hibited explicitly. ^4(^1, ?2, ?3, ^4) is the joint distribution of two quarks and two antiquarks. 
If we are to study the two-pion distribution in a jet, then neither of the two pions should be 
the hadronization of thermal pions only. That is, each pion should contain at least one shower 
parton in order to qualify as a part of the jet. Using the terminology "thermal hadrons" to 
refer to the hadronization of thermal partons only, then in the experimental analysis of jet 
structure such thermal hadrons are regarded as background and are subtracted from the set 
of particles associated with a trigger. In our calculation of dN^-^^^/dpidp2 associated with 
a jet, we simply leave out TT contribution to any pion. Thus there are only two types of 
terms for F4, which we represent schematically as 

F,^{TS){TS) + {TS){SS) , (12) 

where a term of the type {SS) {SS) is omitted because without T it is negligible compared 
to {TS){SS) for hadron pt < 6 GeV/c. In Eq. (12) the parentheses enclose the partons that 
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are to recombine. With that notation {TS){TS) is very different from {TT){SS), which 
contributes to a thermal hadron that we exclude from our consideration. 

In addition to two-pion correlation we shall also study the correlation between a pion 
and a proton in a jet. The single-baryon distribution has the general form 

^odNjB = f (jj F3{qi, q2, gs) Rsiqi, ?2, ?3,p) (13) 
dp J Vi=i Qi J 

where the RF for proton is [8] 

Rp{qi,q2,q3,p) = ^(z/, y2)"+^Z/f+'5 (j^y, - 1 j , = ^ (14) 

g^f) = [B{a + l,a + P + 2)B{a + 1, P + 1)]"' (15) 
with a — 1.75 and (3 — 1.05 [19]. Thus the Trp joint distribution in a high-pr jet is 



1 f f ^ dqi\ 

/ n — ^5(gi, ?2, qs, ?4, 95) Rniqi, q2,Pi)Rp{q3, ?4, ?5,P2)- (16) 

iP2 \i=i qi J 



P\P2dp\dp2 P1P2 
The 5-parton distribution F5 has the schematic form 

F5 = {TS){Trs) + {rs){rss), (17) 

where we have omitted terms of the type {TS){SSS), {SS){TSS) and {SS){SSS) because 
they are all negligible compared to the ones in Eq. (17). 

In the next section we examine in detail the multi-parton distributions F4 and F5 and 
how they contribute to the associated particle distributions, when either the pion or the 
proton is used as a trigger. 
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3 Distributions of Associated Particles 

To select the appropriate two-particle distributions to calculate, let us examine the type of 
quantities that have been measured experimentally. There exist data from RHIC experiments 
that give some properties of charged particles associated with triggers detected in certain 
Pt ranges [6, 7, 14] We shall therefore calculate distributions with similar kinematic ranges. 
Because of detection limitations there is at present no particle identification in either the 
trigger or the associated particles. We shall, nevertheless, do the calculations for specific 
particle species in anticipation for the corresponding data that will become available in 
the future. In particular, we shall consider tt"*" and p as trigger particles, and tt^ and 
p as associated particles. Considering all charged particles would involve overwhelming 
complications and uncertainties without gaining clarity. 

The two-particle distributions given in Eqs. (11) and (16) do not explicitly specify what 
the two particles at pi and p2 are within the same jet. The jet momentum A;, which is 
the momentum of the initiating hard parton, is imbedded in the expression for the shower 
partons, Eq. (2), and is integrated over all values in a collision. When the two-parton 
distribution, such as that given in Eq. (10), is generalized to F4 in Eq. (12), it is important 
to make sure that the hard partons in each of S in (12) are one and the same, and similarly 
in Eq. (17). To make this point explicit, it is helpful to write the two-particle distribution 
within one jet as dNj^^f^^{k) / dpidp2, so that all shower partons are associated with the hard 
parton i at momentum k. Then the particle (/i2) distribution associated to a trigger particle 
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(hi) at pi is 



dNh, 



dpidp2 



(k) 



(18) 



dp2 




{k) 



Further integration of Eq. (18) over pi in the specified range of the trigger momentum yields 
the distribution of the associated particles that is measured. Note that the energy loss 
suppression factor ^ is cancelled in Eq. (18). In the following our expressions for S will not 
contain ^J2il dkkfi{k) that appears in Eqs. (1) and (2), since it is now shown explicitly in 



Let us consider first tt"*" trigger, and tt"*", 7r~, and p associated with it. The initiating hard 
parton i can be u, d, s, u, d, s and g. For every given i, the shower parton j can be u, d, u, 
d. There are therefore 28 possible SPD's. For tt+tt^ in a jet, F4 for four partons has three 
terms 



where the first pair of parentheses in each term correspond to the trigger, the second the 
associated particle. We omit the term {SS){SS) because it is negligible. Note that it is the 
only term contributing to two particles in a jet produced in pp collisions, since there are no 
thermal partons of any significance in the environment. Herein lies already the difference 
between jets produced in heavy-ion and hadronic collisions, without even any quantitative 
details to be investigated. For 7r+7r+, j should only be w or J within each of the six sets of 
parentheses; the thermal partons should just be given the flavor of the complement to make 
a 7r+, i.e., ud. Thus {TS){TS) is a symbohc short-hand for 



Eq. (18). 




■+ 



{rs){rs) + {rs){ss) + {ss){rs) , 



(19) 



{rS){TS) = [{%Ss) + {SuTs)] ■ [{%Ss) + {S^Js)] , 



(20) 
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where the subscript denotes the flavor. For the other terms in Eq. (19) we have 



(21) 



{SS){TS) = (S^S^) [(T„5,-) + {S^T^] . (22) 

For 7r"'"7r~ production we have in a similar way 

Ff"' = [{%Ss) + {S^Ta)] ■ [{%Sd) + {SuTa) + {SuSa)] + (5„5j) [(r„5j) + {S^%)] . (23) 

It should be noted that in the above expression u (d) have momenta (52), and u (d) have 
momenta {q^). 

In the case of vr+p correlation there is an extra T or <S, as shown in Eq. (17). Hence, the 
5-parton distribution has the form 

Ft" = [{%Ss) + iSM] ■ [{TuTuSa) + 2{TA%) + 2{%S^S,) + (<S„<S„T,)] . (24) 

The momenta of the quarks in the proton are in the order u{q3)u{q4)d{q5). The two factors 
of 2 in the above equation arise from the fact that Tu{(lz)Su{Q.A)^d{(l5) makes the same contri- 
bution as Su{q3)%i{q4)Td{q5) . That is, the two u quarks can each receive contributions from 
the thermal and shower sources. The same holds for T^SuSd and SuXiSd also. 

For proton trigger we shall consider only tt"*" and 7r~ as the associated particles. For tt"^ 
associated with p, Fg"^^ is trivially related to simply by interchanging the positions of 
the parentheses in Eq. (24), i.e., interchanging trigger and associated particle, and pi <-> p2- 
For 7r~ associated with p we have 

Fr~ = [{%%Sd) + 2{%Su%) + 2{%SuSd) + {SuSu%)] ■ [{%Sd) + {S^d)] ■ (25) 
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The thermal partons are flavor independent, so Eq. (5) can be used for u, d and their anti- 
quarks. The situation with the shower partons are far more complicated. Their distributions 
depend on the species of the initial hard parton i and the shower parton j. However, by not 
considering the production of kaons and hyperons, there are only three basic SPD's: K, L, 

and G. We can write Sj in the matrix form 

/ \ 
K L 



SI 



L K 
L L 



\ 



i^u,d, s, g, j = u, d. 



(26) 



G G 

For example, 5"^ — K, S'^ — L, Sg — L, Sg — G. Antiquarks are hke quarks of different 
flavors, i.e, = K, = L, = L, etc. Note that hard s quark is included, though not 
shower s quark. The parametrizations of the SPD's are given in Ref. [15], in which one can 
also find a thorough discussion of why j = g is excluded due to gluon conversion. In Eqs. 
(20)-(25) the label i for hard partons is suppressed, but is shown explicitly in Eq. (18). The 
subscripts of S in those equations correspond to the label j above. Thus for every fixed i all 
the Sj distributions in Eqs. (20)- (25) can be rewritten as K, L and G. When i is changed, 
the translation also changes. Thus, for example, when i — u, u, d, d, s, s, and g, the (TS) 
for TT"*" in Eq. (20) becomes 



{rS)i=^ = (TL) + {KT) , 



(27) 



i'rS)i=u,d,s,s - {TL) + (LT) , 



(28) 



(T5),=j = {TK) + (LT) , 
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(29) 



{TS),=, = {TG) + (GT) . (30) 

The index of T has been omitted because of flavor independence. In the case of proton, 
\i i — then only one shower u quark can be valence, the other must be in the sea; 
consequently, we have for the TSS part in Eq. (24), for example, 

{TSS)i=^ = 2{rKL) + {KLT) , (31) 

{rSS),=a = ^{rLK) + (LLT) , (32) 

while, for i — u,d,s and s, K^s above are replaced by L's, and for i — g, both K and L are 
replaced by G. 

The application of Eqs. (27)-(30) to (20)-(22) must take into account the consideration 
that there can only be one valence quark in a jet, and it can be in either the trigger or the 
associated particle. Thus, for example, {TS){TS) in Eq. (20) for tt+tt"^ should have the 
explicit form, ior i — 

[{TS){TS)l^^ = \{[{rL) + {Kr)].[{TL) + {LT)] 

+ [(TL) + (LT)] . [(TL) + {KT)]} , (33) 

whereas, for 7r+p, the term {TS){TSS) in Eq. (24) becomes 

[{rS){rSS)l^^ = \ {[(TL) + {KT)] . [2{TLL) + (LLT)] 

+ [(TL) + (LT)] ■ [2{TKL) + {KLT)]} . (34) 

So far the complications above are due to the differences in SPD's for different i and 
j. Further complications arise when the constraint due to momentum conservation is to 
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be applied. Since the sum of the momenta of shower partons in a jet cannot exceed the 



hard parton momentum k, the momentum fractions in the arguments of the SPD's cannot 
be independent. In Eqs. (12) and (17) there are terms where S appears twice or thrice. 
Let the momenta of the two-quark case be denoted by and q^, and let us consider the 
term (KT) + (LT) in Eq. (33) for illustration. The constraint applies only to K{za) and 
L{zb), where Za^ — Qa,b/k, while the momenta of T are independent. If K has the leading 
momentum g^, then the maximum momentum oi L is k — Qa, and vice- versa. Thus we use 
the symmetrized combination given in Eq. (3) 



{K{Za),L{Zb)}^- 



(35) 



In the case when there are three SPD's, as in Eq. (34), we symmetrize as follows 



{K{Za),L{Zb),L{z,)} = 



K(za)ai-^],L^ 



1 - Zn 



1 - Za 



Zb 



(36) 



The procedure for our calculation is now completely specified. 



4 Results 

To provide an understanding of the order of magnitude of the various terms, let us first 
give the result of tt'^tt''" correlation in central Au-|-Au coUisions at ^/s — 200 GeV. Since the 
STAR data [4] are for the trigger momentum in the range 4 < < 6 GeV/c, we calculate 



dp2 h ^ dp2 



(37) 

7r+(pi) 
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where the RHS is defined by Eq. (18). The result is shown in Fig. 1. The dashed fine indicates 
the contribution from the first term in Eq. (19) that has the structure {TS){TS), while the 
dash-dot fine represents the next two terms in Eq. (19) that are of the form {TS){SS) + 
{SS){TS). The solid line is their sum. Clearly, the overall distribution is dominated by the 
component that involves two thermal partons due to the high density of those soft partons. 
Each time a thermal parton is replaced by a shower parton in their recombination, the yield 
is lower. For that reason we have not bothered to calculate the contribution from {SS){SS), 
which corresponds to the double application of the fragmentation function. Presumably such 
contributions can become important at very large px, where the effects of thermal parton 
are insignificant. To generate four shower partons in a jet resulting in two pions each with 
Pt > 4 GeV/c would require a much harder collision than is necessary if some thermal 
partons can participate. The issue here is not which channels a hard parton hadronizes into, 
given a value of k (as one considers in fragmentation), but rather, given two pions at pi and 
P2, what the most favorable value of k is in the environment of dense thermal partons (as 
one considers in recombination). 

We next consider the dependence on the density of thermal partons. Since we have 
already investigated d+Au collisions in connection with the Cronin effect [12], where the 
soft parton distributions (called thermal also) have been determined for various centralities, 
we calculate tt'^tt"^ correlation for three cases: central Au+Au (0-5%), central d-l-Au (0- 
20%) and peripheral d-l-Au (60-90%). The result is shown in Fig. 2, where the three cases 
are represented, respectively, by solid, dashed and dash-dot fines. Evidently, the density of 
thermal partons has a crucial effect on the yield of the associated particles. Since peripheral 
(60-90%) d-|-Au collisions arc almost equivalent to pp collisions, we can see directly from Fig. 
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2 that the structures of jets in central Au+Au and pp coUisions are drastically different. If 
we plot the ratio of the spectra for central Au+Au to peripheral d+Au in linear scale, we 
get the solid curve in Fig. 3. The ratio of central d+Au to peripheral d+Au is shown by the 
dashed line. The former ratio exceeds 5 around px — 2 GeV/c and remains large throughout 
the intermediate pt range. We now have strong evidence that the structure of jets in nuclear 
collisions is very different from that in hadronic collisions. 

We now consider tt"*", tt" and p production associated with a tt"*" trigger in the 4 to 6 
GeV/c range. Fig. 4 shows the three contributions by the thin solid (7r+), dashed (tt") and 
dash-dot lines (p) , with the sum indicated by the thick solid line. The data points are from 
STAR, which includes all charged hadrons in both the trigger and the associated particles 
[6]. Since what is calculated is not exactly what is measured, one should not expect perfect 
agreement. However, tt'^ is a dominant component of the trigger, and the data average 
over different trigger particles, whereas the different associated particles are summed. Thus 
what is calculated should not differ greatly from what is measured. Indeed, the agreement 
is very good both in normalization and in shape. It is therefore reasonable to infer that our 
approach has captured the essence of the physics of hadronization. 

It is interesting to note that the yield of 7r~ is higher than that of tt"*" when the trigger 
is TT"*", as is evident in Fig. 4. The reason is that when i is summed over all hard parton 
species, the cases when i — d and u can give rise to valence shower partons that enhance the 
7T~ production through the K distribution, but not the tt"*" production. When i = u and d, 
the trigger uses up the valence shower parton to form tt"*" so the associated particle, whether 
TT"^ or 7r~, has to be formed by the sea shower parton through the L distribution, resulting in 
no big difference between 7r+ or . Thus adding up the contributions from all hard partons 
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results in more 7r~ than tt^ in a tt"^ triggered jet. This is a prediction that can be checked 
by experiments with good particle identification. 

We also note that the proton yield in Fig. 4 is greater than the n'^ yield in the 2-3 GeV/c 
range because a u or d shower parton can recombine with two thermal partons, thereby 
increasing the p/ir ratio for the same reason that the ratio exceeds 1 without trigger [8]- [11]. 
However, the proton yield is less than the Tr~ yield in the tt'^ triggered jet, since ir'^p does 
not have the advantage of the zZ-initiated jet that enhances the tt+tt" production. 

Finally, we come to the proton trigger and show in Fig. 5 the result of our calculation 
for the associated particles being 7r+ and 7r~, in dashed and dash-dot lines, respectively. 
The solid line is their sum. The normalization and shape of the total distribution for the 
associated particles are roughly the same as those in Fig. 4 for 7r+ trigger. The 7r+ and tt" 
components have no noticeable difference (bearing in mind that the 7r~ curve is lowered by a 
factor of 2 to avoid overlap). That is reasonable, since i — d and u favor tt"*" and 7r~ equally, 
while i — u and d are both used in the trigger, leaving tt"^ and 7r~ again on comparable footing. 
These features can also be checked directly by experiments when particle identification is 
improved. 

In the foregoing we have presented the distributions of the associated particles, which 
offer more details than the overall yields. The latter provide a short and useful summary of 
the jet structure that is easier to measure. We therefore calculate the three lowest moments 
of the distributions that have already been obtained 



where dN^^^ydp2 is the pr distribution of all the particles associated with trigger hi that 




(38) 
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we have calculated: tt"*", tt" and p for hi — n'^, and tt"^ and 7r~ for hi — p. The lower limit 
of the integral in Eq. (38) is set at 0.5 GeV/c because our calculated result is not reliable 
for pt < 0.5 GeV/c; the upper hmit is set at 4.5 GeV/c, since we do not want it to exceed 
the average of the trigger momentum that is between 4 and 6 GeV/c. Thus, by definition, 
is a measure of the average number of particles associated with trigger hi, Mi^^^ being 
the total scalar pt of those particles, and M^^^'' the total p^ of them. The last quantity 
is insensitive to the low p2 behavior of dN^'^^y dp2, and is a good measure for comparison 
between theory and experiment. Our results on M^'^^^ for central Au+Au coUisions are shown 
in Table I. Note that the values of the three moments change by mildly decreasing factors 

Table 1: Values of the moments M^^^^ for central Au+Au at ^/s = 200 GeV 



n 





1 


2 




1.394 


1.707 


2.703 


hi=p 


0.882 


0.999 


1.450 



(0.63, 0.59, 0.54) when the trigger is changed from tt"^ to p. This is a feature that can more 
easily be checked by experiments than the distributions dN^'^^ydp2 themselves. 

To compare the above results with what one can expect from pp collisions, we can return 
to what we have already calculated, i.e., tt+tt"*" correlation in peripheral d+Au collisions, since 
that is very close to pp collisions. Indeed, for an appreciation of the centrality dependence 
we calculate the moments of the three distributions shown in Fig. 2. The results are given 
in Table II. The ratios R = [Au+Au (central)]/[d+Au (peripheral)] are 4.8, 4.9 and 4.9 for 
n = 0, 1, 2, respectively. They are very similar and roughly 5. We expect that if all particles 
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Table 2: Values of the moments M^^^^^ for central tt"*" associated particle only in three 
colliding systems 



n 





1 


2 


Au+Au (central) 


0.428 


0.498 


0.754 


d+Au (central) 


0.200 


0.223 


0.331 


d+Au (peripheral) 


0.089 


0.101 


0.153 



associated with the trigger are included, the ratio will remain about the same. 

There are some data on total charged multiplicity and total scalar p^, but they include 
the trigger [6]. The ratio for [Au+Au (central)]/]?]? on multiplicity is ~ 1.3, and on scalar pr 
is ~ 1.5. By subtracting out the trigger contribution, it is possible to see a rough agreement 
with our result; however, we leave the quantification of the comparison to the experiments. 

5 Conclusion 

We have studied the structure of jets produced in heavy-in coUisions by calculating dihadron 
correlation in the framework of parton recombination. Since the jets are produced in the 
environment of dense partonic medium, they are different from the ones produced in pp 
collisions and e^e" annihilation. The interaction between the hard and soft partons is very 
important. It is taken into account in our study by allowing shower partons to recombine 
with the thermal partons. Since our formalism has been applied successfully to single-particle 
spectra in previous studies, we have no freedom to adjust any part of our treatment of the 



18 



two-particle distributions, nor is there any free parameter to vary. All the results shown in 
the previous section are predictions. 

In our approach to hadron production at high and intermediate pr the effect of energy loss 
by the hard partons traversing the dense medium is represented by a multiplicative factor ^. 
That factor is cancelled in our definition of the associated particle distribution, which is the 
ratio of the two-particle distribution to the trigger-particle distribution. Thus the dihadron 
correlation in a jet that we calculate is independent of the degree of jet quenching. It can be 
compared to the corresponding data on the near-side jet, since such jets are produced by hard 
collisions near the surface facing the detector and suffer minimal energy loss. Unfortunately, 
inadequate particle identification renders unfeasible direct comparison between the currently 
available data with our predictions. 

In principle, it is possible to calculate what has currently been measured, i.e., all charged 
hadrons in the trigger and other particles in the jet. In practice, the task would be dauntingly 
complicated and involve many more terms than what we have already considered, including 
thermal s partons whose parametrizations have been less reliably determined. Besides, the 
current data are in a passing phase; better particle identification is forth-coming. What 
we have calculated are for clean triggers and associated particles, and can be effectively 
compared with future data. 

It is evident from the results of our study that the dihadron correlation is dominated 
by the components that involve the highest number of thermal partons. The recombination 
mechanism boosts the yield when high-density thermal partons are included, but also boosts 
the pt of the product when the semi-hard shower partons are involved. The effect has 
already been shown to operate in the singlc-particlc distributions, but now exhibits itself 
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more conspicuously in dihadron correlations in jets, since at least four partons are involved, 
two of which can be thermal. The difference between jets produced in heavy-ion collisions 
compared to those produced in pp collisions is huge, as evidenced by the large ratio shown 
by the sohd hne in Fig. 3 and by the ratios of the moments in Table II - about 5 between 
the first and third rows of values. 

It is hard to see how the properties of the dihadron correlations in jets that we have 
found can be reproduced in any fragmentation model even with medium modification of the 
fragmentation function, whose focus has been on the effect of energy loss [20]. There are 
terms in our recombination formulas, when the shower parton momenta are symmetrized, 
that cannot be written in factorizable forms involving products of two FF's. Leaving aside 
such technical details, let us accept the conceivable possibility that perturbative branching 
of a hard parton can generate hard shower partons, and that they can further fragment by 
means of suitable modified FF's that mimic TS recombination. If it is successful in the end, it 
seems that the scheme would have lost the original advantage of the fragmentation approach 
that relies on the universality of the FF's. If the modification of a FF depends sensitively 
on the detailed properties of the partonic medium, then, by comparison, TS recombination 
would seem to be a more direct and physically cogent approach to hadronization. Besides, 
as pointed out earlier, energy loss is not the issue in dihadron correlation. 

There are hmitations to the formahsm that we have used for our calculations. We have not 
considered the dependence of dihadron correlation, since the SPD's used are for Q fixed 
at 10 GeV/c [15]. That hmitation is not a matter of principle, but of practice. To account 
for the evolution of the SPD's is a worthwhile problem in its own right, inasmuch as the 
same problem for the FF's has been pursued for decades [21, 22]. To apply that evolutionary 
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property to the dihadron correlation would be prohibitively complicated. Attempts have been 
initiated to investigate the evolutionary aspect of dihadron distribution of the fragmentation 
process in the operator formalism [23]. An input on the initial distribution for such an 
evolution would still have to involve the type of consideration presented here. Furthermore, 
how the thermal partons are to be incorporated in that approach is not clear. 

Another limitation is rooted in our formalism. Since we have relied entirely on our ID for- 
mulation of recombination, it is not possible in the same formalism to address the question 
of angular correlation in a jet. Since jets are 3D objects, it is obvious that longitudinal- 
transverse correlation within a jet can contain information about the properties of the re- 
combining partons that we cannot probe in the simplified ID formalism. Clearly, despite 
the substantial progress that has been made in our line of investigation, there remains much 
ground to improve and generalize in the study of the physics of hadronization in heavy-ion 
collisions. 

Acknowledgment 

We are grateful to X. N. Wang, A. Majumder, F. Wang, A. Sickles and M. Tannenbaum 
for helpful discussions and communication. This work was supported, in part, by the U. 
S. Department of Energy under Grant No. DE-FG02-96ER40972 and by the Ministry of 
Education of China under Grant No. 03113. 



21 



References 

[1] K. Adcox et al, (PHENIX), Phys. Rev. Lett. 88, 022301 (2002); C. Adler et al, (STAR), 
Phys. Rev. Lett. 89, 202301 (2002); S.S. Adler et al, (PHENIX), Phys. Rev. Lett. 91, 
072301 (2003); J. Adams et ai, (STAR), Phys. Rev. Lett. 91, 172302 (2003). 

[2] D. d'Enterria, talk given at Quark Matter 2004, nucl-ex/0404018. 

[3] For a recent review see P. Jacobs and X. N. Wang, hep-ph/0405125. 

[4] C. Adler et al, (STAR), Phys. Rev. Lett. 90, 082302 (2003); J. Adams et al, (STAR), 
Phys. Rev. Lett. 91, 072304 (2003). 

[5] X. N. Wang and M. Gyulassy, Phys. Rev. Lett. 68, 1480 (1992); X. N. Wang, Phys. 
Rev. C 58, 2321 (1998); 61, 064910 (2000). 

[6] F. Wang, (STAR), nucl-ex/0404010. 

[7] J. Rak (PHENIX), nucl-ex/0306031. 

[8] R. C. Hwa, and C. B. Yang, Phys. Rev. C 67, 034902 (2003). 

[9] V. Greco, C. M. Ko, and P. Levai, Phys. Rev. Lett. 90, 202302 (2003); Phys. Rev. C 
68, 034904 (2003). 

[10] R. J. Fries, B. Miiller, C. Nonaka and S. A. Bass, Phys. Rev. Lett. 90, 202303 (2003); 
Phys. Rev. C 68, 044902 (2003). 

[11] R. C. Hwa and C. B. Yang, nucl-th/0401001, Phys. Rev. C (to be pubhshed). 

[12] R. C. Hwa and C. B. Yang, nucl-th/0403001, Phys. Rev. Lett, (to be published) 

22 



[13] R. C. Hwa and C. B. Yang, nucl-th/0404066, Phys. Rev. C (to be published). 
[14] A. Sickles, (PHENIX), nucl-ex/0403028. 

[15] R. C. Hwa and C. B. Yang, hep-ph/0312271, Phys. Rev. C (to be published). 

[16] K. P. Das and R. C. Hwa, Phys. Lett. 68B, 459 (1977). 

[17] R. C. Hwa, Phys. Rev. D 22, 1593 (1980). 

[18] R. C. Hwa and C. B. Yang, Phys. Rev. C 66, 025205 (2002). 

[19] R. C. Hwa, and C. B. Yang, Phys. Rev. C 66, 025204 (2002). 

[20] X. N. Wang, Z. Huang and I. Sarcevic, Phys. Rev. Lett. 77, 231 (1996); X. N. Wang 
and Z. Huang , Phys. Rev. C 55, 3047 (1997). 

[21] J. Binnewies, B. A. Kniehl and G. Kramer, Z. Phys. C 65, 471 (1995); B. A. Kniehl, 
G. Kramer and B. Potter, Nucl. Phys. B597, 337 (2001). 

[22] S. Kretzer, Phys. Rev. D 62, 054001 (2000). 

[23] A. Majumder and X. N. Wang, hep-ph/0402245. 



23 



Figure Captions 

Fig. 1. Transverse momentum distribution of tt"'" associated with a tt^ trigger. The contri- 
bution from terms of the form {TS){TS) is shown in dashed hne, while the contribution 
from terms of the form {TS) {SS) is shown in dashed-dot hne. The sohd hne is the 
sum of aU components. 

Fig. 2. Associated particle (tt"*") distribution with 7r+ trigger for central Au+Au (solid), 
central d+Au (dashed) and peripheral d+Au coUisions (dash-dot line). 

Fig. 3. The ratio of central Au-|-Au to peripheral d-l-Au coUisions (solid line) and that of 
central d-l-Au to peripheral d-l-Au collisions (dashed line). 

Fig. 4. Transverse momentum distributions of tt"*", tt" and p, associated with a 7r+ trigger. 
The data are from STAR [6] for all charged hadrons in the trigger and associated 
particles. 

Fig. 5. Transverse momentum distributions of tt"*" and 7r~ (lowered by a factor of 2) asso- 
ciated with proton trigger. 
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